In vivo anomalous diffusion and weak ergodicity breaking of lipid granules 



O 

(N 



I 

C 

O 
o 



(N 
> 

o 



X 



Jae-Hyung Jeon, 1 Vincent Tejedor, 2 Stas Burov, 3 Eli Barkai, 3 Christine 

Selhuber-Unkel, 4, 5 Kirstine Berg-S0rensen, 6 Lene Oddershede, 4 and Ralf Metzler 1 

Physics Department T30g, Technical University of Munich, 85747 Garching, Germany 
2 Physique Theorique de la matiere condensee, Universite Pierre et Marie Curie, 4 place Jussieu, 75252 Paris, France 
3 Physics Department, Bar-Ilan University, Ramat Gan 52900, Israel 
^ Niels Bohr Institute, Blegdamsvej 17, 2100 K0benhavn 0, Denmark 
5 Institute for Materials Science, University of Kiel, Kaiserstrafie 2, 24143 Kiel, Germany 
6 Physics Department, Technical University of Denmark, 2800 Kongens Lyngby, Denmark 

(Dated: 18th January 2011) 

Combining extensive single particle tracking microscopy data of endogenous lipid granules in 
living fission yeast cells with analytical results we show evidence for anomalous diffusion and weak 
ergodicity breaking. Namely we demonstrate that at short times the granules perform subdiffusion 
according to the laws of continuous time random walk theory. The associated violation of ergodicity 
leads to a characteristic turnover between two scaling regimes of the time averaged mean squared 
displacement. At longer times the granule motion is consistent with fractional Brownian motion. 
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Crowded colloidal systems of hard and soft core 
particles display a rich physical behavior including dy- 
namic arrest, non-exponential relaxation, and anomal- 
ous diffusion In biological cells the cytoplasm is a 
crowded soft core colloid made up of large biopolymers 
such as ribosomes, proteins, or RNA, occupying volume 
fractions of 34% or above @, Q , amidst small particles 
such as water, ions, and lower mass fraction biopolymers. 
Such crowding significantly impacts cellular biochemical 
reactions 0] and effects subdiffusion of the form 0, [B| 



(r 2 (t)) (xK a t a , 0<a<l, 



(1) 



of larger molecules or tracers in living cells 
While normal diffusion, by virtue of the central limit the- 
orem, is characterized by the universal Gaussian probab- 
ility density function and therefore uniquely determined 
by the first and second moments [l7[ , anomalous diffusion 
of the form ([1]) is non-universal and may be caused by 
different stochastic mechanisms. These would give rise to 
vastly different behavior for diffusional mixing, diffusion- 
limited reactions, signaling, or regulatory processes. To 
better understand cellular dynamics, knowledge of the 
underlying stochastic mechanism is thus imperative. 

Here we report experimental evidence from extensive 
single trajectory time series of lipid granule motion in 
Schizosaccharomyces pombe (S. pombe) fission yeast cells 
obtained from tracking with optical tweezers (resolving 
10~ 6 sec to 1 sec) and video microscopy (10 -2 sec to 
100 sec). Using complementary analysis tools we demon- 
strate that at short times the data are described best 
by continuous time random walk (CTRW) subdiffusion, 
revealing pronounced features of weak ergodicity break- 
ing in the time averaged mean squared displacement. At 
longer times the stochastic mechanism is closest to sub- 
diffusive fractional Brownian motion (FBM). The time 
scales over which this anomalous behavior persists is rel- 



evant for biological processes occurring in the cell. An- 
omalous diffusion may indeed be a good strategy for cel- 
lular signaling and reactions 0, HJ . 

CTRW and FBM both effect anomalous diffusion of 
the type ([T]) [18| . Subdiffusive CTRWs are random walks 
with finite variance (Sx 2 ) of jump lengths, while the 
waiting times between successive jumps are drawn from 
a densit y ib (t) ~ T a /t 1+a with diverging characteristic 
time [l,ll3|- Such scale- free behavior results from mul- 
tiple trapping events in, e.g., comb-like structures or 
random energy landscapes [2fJ . Power-law waiting time 
distributions were also identified for tracer motion in re- 
constituted actin networks |2l|. Subdiffusive FBM is a 



random process driven by Gaussian noise £ with long- 
range correlations, (£(0)£(t)) ~ a{a - l)t a ~ 2 [H], and it 
is related to fractional Langevin equations [23| . The finite 
characteristic time scales associated with FBM contrast 
the ageing property of the subdiffusive CTRW processes 
[Supplementary Material (SM)]. 

Single particle tracking microscopy has become a 
standard tool to probe the motion of individual tracers, 
especially inside living cells, and provides valuable in- 
sights into cellular dynamics [6j . The recorded time series 
r(i) are analyzed in terms of the time averaged mean 
squared displacement as function of the lag time A, 



S 2 (A,T) 



1 



T — A 



T-A 



r(t + A)-r(t) dt, (2) 



where T is the total measurement time. CTRW and 
FBM subdiffusion show markedly different behavior of 
S 2 (A,T) [US]. Subdiffusion of the form S 2 ~ A" 
was found in several tracking experiments in vivo 0- 
Hij ]. Crowding- induced subdiffusion was also observed 
in controlled environments , e.g., in concentrated protein 
and dextran solutions 16, 28|-30|]. While FBM was pro- 
posed as stochastic mechanism in some of these systems 
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Figure 1: Time averaged mean squared displacement 8 2 from 
individual trajectories of lipid granules in S.pombe in Early 
Mitotic cells (lower curves) and in Early Telophase (up- 
per curves), measured by optical tweezers 33]. A distinct 
turnover from P ~ A to ~ A' 3 (j3 ns 0.10... 0.20) occurs. 
The two thick lines in the middle show the averages of the 
ET (red O) an d EM (black □) data sets. The overlaid thick 
black line is the result of CTRW simulations in an harmonic 
potential. ET versus EM curves are shifted vertically. 



14 Il6l l31j CTRW subdiffusion has not yet been identi- 
fied. Here we demonstrate that the short time motion of 
lipid granules in S.pombe cells follows the laws of CTRW 
subdiffusion and features weak ergodicity breaking. 

Short time behavior. Granule motion was recorded in 
an optical tweezers setup (SM). In the experiment the 
trap is initially centered onto the granule such that no 
force is exerted. When the granule starts to move away 
from the trap center it experiences a restoring Hookean 
force |32j . The measured S 2 (A) curves from two different 
cell stages are shown in Fig. [T] [33| , along with a typical 
sample trajectory in Fig. 5a (SM). In Fig. [T] a distinct 
turnover is observed from an initial linear growth S 2 ~ A 
to a power-law behavior S 2 ~ A^ with f3 ~ 0.15 . . . 0.20. 

The size of the lipid granules is about 300 nm. In addi- 
tion to the amplitude scatter between different trajector- 
ies expected from CTRW theory 3, 2(| , the fluctuations 
of the data in Fig. [T]are due to natural granule size vari- 
ations and different optical conditions for each trajectory. 
We note that oscillations around the turnover may arise 
for subdiffusion in an underdamped medium [34j]. 

Evaluating extensive data for the granule motion ana- 
logous to the data presented here by ensemble averages, 
anomalous diffusion with a « 0.80... 0.85 yields in a 
range of 0.1 to 3 msec 13] . In contrast, in Fig. [1] the 



initial behavior of S 2 corresponding to this time range 
does not exhibit any apparent anomaly but scales like 
S 2 ~ A 1 . At longer lag time A due to the trap force one 
would expect 5 2 to saturate to a stationary thermal value; 
instead, the regime S 2 ~ A* 3 appears. This behavior is 
consistently observed in different cell stages (Fig. [T|). 



Such a peculiar behavior is fully consistent with CTRW 
subdiffusion: For free motion one finds S 2 ~ A/T 1- ", 
whose A scaling is independent of a, while the corres- 
ponding ensemble average follows Eq. ([1]) [iH, [25| . Under 
confinement a turnover to the power-law S 2 ~ (A/T) 1- " 
occurs [HI23. This second power-law regime is termin- 
ated when A approaches the total measurement time T, 
causing a dip in S 2 back to the plateau of the ensemble 
average. The observed characteristic turnover behavior 
is intimately connected to CTRW ageing and ergodicity 



breaking [35J, |36| . Corralled motion [37[ could not explain 
the observed behahior 5 2 ~ A turning over to S 2 ~ A^. 

The following features unanimously point toward 
CTRW subdiffusion as the stochastic mechanism for the 
granule motion at short times: (i) The time average S 2 
initially scales linearly with A, albeit the ensemble aver- 
age shows subdiffusion of the type ([!} in comparable time 
ranges, (ii) At longer times a turnover to the power-law 
S 2 ~ A' 3 occurs instead of the convergence to a plateau, 
which would necessarily occur for an ergodic process; the 
anomalous diffusion exponent a rj 0.80 . . . 0.85 observed 
in the ensemble average is consistent with the exponents 
f3 w 0.15 . . . 0.20 observed in Fig. Q] based on the relation 
(3 = 1 — a [2f| [I?} , as well as the slopes of the long time 
data (see below), (iii) The specific form of 5 2 of the data 
nicely coincides with simulations results for a subdiffusivc 
CTRW in an harmonic potential (Fig. [I). These observa- 
tions prove that the granule motion in fact exhibits weak 
ergodicity breaking associated with subdiffusivc CTRW 
processes. This is the central finding of this work. 

We rule out the possibility that the initial CTRW 
simply turns over to much slower diffusion with a ~ 0.15 
. . . 0.20, as the overlapping video tracking data shows 
subdiffusion with a s=s 0.8 or slightly below. Also, CTRW 
or FBM trajectories with a = 0.20 would cause a quite 
extreme stalling or antipersistence, respectively, which is 
inconsistent with the recorded time series [Fig. 5 (SM)]. 

What results yield from complimentary criteria? (i) 
Evaluating <5 2 (A, T) as function of the total measurement 
time T, no ageing is observed [Fig. 6 (SM)], contrasting 
the scaling 6 2 {T) ~ T a_1 predicted for CTRW subdif- 
fusion with diverging mean waiting time (t) [13-0] . (ii) 
Determining the distribution </>(£) of the deviations of 
the relative time averaged mean squared displacement 
£ = 5 2 1 (S 2 ) around the ensemble mean (£ = 1) we find 
the bell-shaped form of Fig. [2] with (f)(0) « 0. This form 
deviates from the predicted shape for CTRW subdiffu- 
sion with (t) = 00 where </>(0) ^ [2^|, see Fig. [5] 

These seemingly conflicting observations can in fact 
be reconciled when we consider a power-law waiting time 
distribution with a cutoff effective at times t t* , 



dt 



1 - 



(t + t) a 



exp 



(3) 



Physically such cutoffs occur naturally in finite systems, 
corresponding to, e.g., a maximal well depth in a random 
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Figure 2: Distribution <^>(£) of the relative deviations £ = 
S 2 /{8 2 ) of J 2 obtained from averaging the data of Fig. [1] (Early 
Telophase, Early Mitosis). Simulations results are shown 
based on Eq. © with r* = 182 msec (□), and for a CTRW 
with (t) = oo (O); with parameters a = 0.85, T = 3.0 sec, 
r = 0.04545 msec, and A = 0.04545 msec. See Fig. 7 (SM). 



Figure 3: Time averaged mean squared displacement for a 
CTRW with waiting time cutoff (|3]) in an harmonic potential. 
Inset: Comparison of the average behavior of the simulated 
trajectories ( — ) and the average of the experimental data 
from Fig. [TJ (Q, D)- Parameters as in Fig. [2] 



energy landscape. In free space Eq. ([3]) produces initial 
subdiffusion (r 2 (i)) ~ t a , turning over to (r 2 (t)) ~ t 
at i > r*. As shown in Fig. |3] the cutoff- waiting time 
distribution (j3|) reproduces well the experimental beha- 
vior from Fig. [TJ The scatter distribution 0(£) [Figs. 
7 (SM)] shows excellent agreement, and the ageing plot 
[Fig. 6 (SM)] is consistent with the data. Here we chose 
the cutoff time t* = 182 msec, and t = 0.04545 msec. 

Collecting all results we conclude that at short times 
the granule motion is described by CTRW subdiffu- 
sion with truncated power-law waiting time distribu- 
tion ([3]). While the cutoff ensures that the typical scal- 
ing of S 2 (A,T) still exhibits the weak ergodicity break- 
ing features, for appropriate choice of the cutoff time 
it is consistent with the absence of ageing effects, i.e., 
5 2 {T) ~ T°, and the observed form for 

At longer times the motion of the lipid granules was 
recorded by video particle tracking. Fig. [4] shows the 
time averaged mean squared displacement. Initially the 
slope is around a w 0.8 or slightly below, consistent with 
the short time data. Several of the curves turn to a 
gentler slope at around 100 msec, some curves eventually 
switch to normal diffusion (a = 1) at 1 sec. Significant 
deviations are observed, which, for a living systems, is 
not surprising. Both the granule size and the materi- 
als properties of the cellular environment may change on 
these time scales (depolymerization/repolymerization of 
the cytoskclcton etc). As shown in Fig. 8 (SM) the dis- 
tribution of (/>(£) is bell-shaped; no ageing is observed. 
Analysis of the moment ratios for both normal moments 
and the mean maximal excursion statistics (see Ref. [IH 
for details) are consistent with FBM for the range of 
a ss 0.80 . . . 0.85 [Fig. 9 (SM)], as are typical antipersist- 



ent trajectories [Fig. 10 (SM)]. Moreover, the scaling ex- 
ponents of the mean maximal excursion second moment 
are consistently above those of the corresponding regular 
second moments, as predicted for FBM [38|. The velocity 
autocorrelation of the data is consistent with both FBM 
and CTRW [Fig. 11 (SM)], and therefore is not conclus- 
ive. The p- variation method [31] was not conclusive. 

Our analysis of extensive single particle tracking data 
of lipid granules in S.pombe cells demonstrates that at 
shorter times in the msec range the motion displays sig- 
nificant effects of weak ergodicity breaking both in the 
limit of free motion and in the presence of the restor- 
ing trap force. Concurrently no ageing occurs and the 
distribution of time averages 0(£) is bell-shaped around 
£ = 1. We showed that these features are consistent 
with a CTRW process with truncated power-law waiting 
time distribution. At longer times the motion is best 
described by subdiffusive FBM, although a conclusive 
statement in this time range is more difficult due to the 
fact that cellular processes appear to be superimposed 
on the motion. Physically, the CTRW-like motion may 
be associated with the interaction between granules and 
the semiflexible filaments of the cytoskeleton similar to 
the observations in Ref. |2l|. While the behavior be- 



comes more erratic in the long time data a turnover to a 
gentler slope is observed before a final increase to normal 
diffusion. This behavior may be connected to the vis- 
coelastic properties of the complex cellular environment. 
The identification of FBM as stochastic mechanism is 
consistent with conclusions in Refs. 
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Lipid granules provide a natural, inert tracer to explore 
the diffusion properties inside living cells. Subdiffusion 
of large biopolymers in the cell supports the emerging, 
more local picture of cellular transport and regulation 
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Figure 4: Time averaged mean squared displacement of gran- 
ules in S.pombe (cells in interphase) from video tracking data 
(SM). An initial slope around a ~ 0.8 is found, turning over 
to a gentler slope. Several trajectories later exhibit a ~ 1.0. 



CTRW gives rise to a dynamic spatial localization 
1351 ] , while FBM provides a more compact spatial ex- 
ploration and thus increases the local encounter probab- 
ility Q. Identifying these different subdifFusion mechan- 
isms in cells can give new insights into in vivo molecular 
processes at different timescales. 
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